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Abstract 



H 
in 

"jrt \ Consider an exponential dispersion model (EDM) generated by a 

probability /x on [0, oo) which is infinitely divisible with an unbounded 
Levy measure v. The Jorgensen set (i.e., the dispersion parameter 
space) is then M"*", in which case the EDM is characterized by two 
^ ' parameters: ^o the natural parameter of the associated natural expo- 

'^ . nential family and the Jorgensen (or dispersion) parameter t. Denote 

*vq \ by EDM{6Q,t) the corresponding distribution and let Yt is a r.v. with 

cn ' distiihution EDM {6o,t). Then if z^((x,oo)) ~ —£ log x around zero we 

l/-^ . prove that the limiting law Fq of Y^^ as t — )■ is of a Pareto type 

^^ I (not depending on ^o) with the form Fq{u) = for u < 1 and 1 — u~ 

for u > 1. Such a result enables an approximation of the distribu- 
tion of Yt for relatively small values of the dispersion parameter of the 
corresponding EDM. Illustrative examples are provided. 
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1 Introduction 

Let {Ft : < t < to < 00} he a family of distributions associated with posi- 
tive r.v.'s {Yt} and Laplace transforms (LT's) Lt{u) = £'(e~"^'). Also let Yq 
be a r.v. with distribution Fq. Bar-Lev and Enis (1987, Theorem 1) showed 

that Fj"* — )■ yj] as t — 7- (where — )■ designates a convergence in distribution) 
iff 

rimLtiu'/') = Foiu) = l-Foiu) (1) 

at all continuity points of Fq. As Bar-Lev and Enis indicated, such a re- 
sult can be viewed as a "centralization" problem in the following sense. In 
many cases the limiting distribution of Y^ as t — )■ is degenerate. A mea- 
surable transformation gtiYt) is then sought whose limiting distribution is 
non-degenerate. Accordingly, their Theorem 1 suggests a consideration of 
gtiXt) = Yf^ (or, equivalently, of —t\n.Yt) whose limiting distribution is 
non-degenerate (provided that ([T]) is satisfied). Bar-Lev and Enis presented 
several examples which satisfy ([1]). However, these examples heavily depend 
on the explicit (and relatively 'nice') form of Lt. 

A natural question then arises: Can one delineate subclasses of distribu- 
tions which satisfy ([1]), regardless of the explicit form of Lf? Indeed, in this 
note we provide such a subclass, namely a subclass of exponential dispersion 
models (EDM's) generated by a probability /x on [0, cxd) which is infinitely 
divisible of type 1 (c.f., Jorgensen, 1987, 1997, 2006, and Letac and Mora, 
1990). For such a subclass, the EDM is characterized by two parameters: 9q 
the natural parameter of the associated natural exponential family (NEF) 
and the Jorgensen (or, equivalently, the dispersion parameter) t G M"*". We 
denote such a subclass of distributions by EDM{Oq, t) and prove that if Yt has 

a distribution in EDM{9o, t) then Fj~* — )■ Yq as t — )■ 0, where the distribution 
Fq of Yq is of a Pareto type (not depending on 6*0) with the form 

, . jo, if M < 1, 

for some £ > 0. Such a result enables an approximation of the distribution of 
Yt for relatively small values of the dispersion parameter of the corresponding 
EDM. 

This note is organized as follows. In Section 2 we first introduce some 
preliminaries on NEF's and EDM's and then present our main result. Section 
3 is devoted to some examples. 



2 Preliminaries and the main result 

We first briefly introduce some preliminaries related to NEF's and their as- 
sociated EDM's. Let // be a probability measure on M. Assume that the 
effective domain of fi has a nonempty interior, i.e., 

e = mt D = loeR: L{d) = f e-^"" ^i{dx) < oo l 7^ (/>. 

The NEF generated by /x is then given by the set of probabilities 

p{e,fx)idx) = —-fi{dx), eeD 

Note that since /i is a probability measure then E D. The Jorgensen set is 
defined by 

A = |t G M^ : L* is a LT of some measure fit on M} , 

whereas the corresponding EDM is the class of probabilities 

P{e,t,fit){dx) = —-fit{dx),ee D,t G A| , (2) 

Note that the class of EDM's is abundant since any probability measure 
with a LT generates an EDM. An EDM is therefore parameterized by the 
two parameters {6, t) G -D x A, where 6 is the natural parameter of the 
corresponding NEF and t is termed the dispersion parameter. EDM's have 
a variety of applications in various areas, in particular, in generalized linear 
models (replacing the normal model as the error model distribution) and 
actuarial studies. Note that if /i is infinitely divisible then the Jorgensen set 
(or, equivalently, the dispersion parameter space) is A = IR+. Also note that 
if Yt is a r.v. with distribution in (|2]) then its LT is given by 

..y„ (U0 + S) 



Now we consider the case where fi is infinitely divisible law of type 1 
concentrated on IR+. By this we mean that there exists an unbounded positive 
measure v on (0, cxd) such that for ^ > one has L{6) = J°° e~^^fi{dx) = e'^^^^ 
with k{6) = — L (1 ^ e~^^)v{dx) and J„ min(l,a;)z/((ia;) < 00. Therefore 



V is the Levy measure of /i. The Levy process associated with such a /x is 
sometimes called a pure jump subordinator (in this respect, of Levy measures 
for NEF's, see also Kokonendji and Khoudar, 2006). Note that this implies 
that lim6i_j.oo k{ff) = — oo since u is unbounded and therefore lim5i_j.oo L{9) = 
and A*({0}) = 0. We are now ready to present our main result relating to the 
limiting distribution of Y^^* as t — ?■ 0. 



Proposition 1 Assume that fi is an infinitely divisible probability measure 
of type 1. Also assume that G{x) = J^ ^{dy) = ^{{x, oo)) is such that 

lim G'(x)/logx = —i 

X— i-O 



for some £ > 0. Let Oq > 0, then 






^5 



implying by ^\) that Fj * — )■ Yq o^s t — > 0, where the distribution Fq of Yq is 
given by 

1 — M , if U > 1, 



Proof. We first prove (^ for ^o = 0. For this, we give another presentation 
of k in terms of G which is obtained by an integration by parts. For e > 
consider the Stieltjes integral 

/oo /"OO 

(1 - e-''")z/(cix) = {e-'^' - l)G{e) -9 e-^''G{x)dx. (5) 

Since for e — >■ we have (e"^" — 1) ~ —9e and G{e) ~ — £loge we get 

k{e) = \imk,{e) = -e / e-^^'G(x)rfx. (6) 

^^0 Jo 

If < M < 1 we have limi^o""^^* = 0. Since /x is a probability L{0) = 1 
and thus linif^o -^(^^^*)* = 1. If ti = 1 we also have that lim(^o-^(l)* = 1- 
If M > 1, We fix an arbitrary < e < i. By the definition of i there exists 



< T] < 1 such that if < x < r^ then —{i — e) logx < G{x) < —{i + e) logx. 
We now use (E]) for writing 



< -eO I e-^^'logxdx (7) 





/•oo rrj 

k{9) + 9 e-^''G{x)dx -W e"''^ log xdx 

Jr] Jo 

and observing that 

/•OO /"OO 

e / e~'^'=G(x)dx < G{r])e / e-^^rfx = G'(r/)e-''^ ^ 0. (8) 

We need now the foUowing evaluation. For r^ > we have 

1 n 

lim - — -6 / 6-^"= log xdx = -1. (9) 

e^c^log^ 7o 

To prove ([9]), we obtain, by a change of variable to w = 6'x, that 

I rv ]^ /•'ys /-^^ 

-^ / e'^"" log xdx = - — - / e'^logvdv — / e~'"dv — )'e-^oo= — 1, 



log Jo log 9 Jo 

(10) 
where the last term on right hand side of ( TTOl) follows since J^ e "" log f (if 
converges and f^ e~^dv = 1. We now divide both sides of (^^ by log^ and 
let 6^ — )■ OO. From ([8]) and ([9]) we get that for all e > 

-i-ei< liminf t^K^) < hmsup ■r^A;(e) < -£ + ei, 
e^oo log 6* e_^oo log 6* 

i.e., limg^oo ^Te^l^) = ~^- Applying this to 6* = m^/* with a fixed -u > 1 and 
letting t — i- 0, we get 

limtfcfn^/*) = -flogM, lim(L(u^/*))* = 4, 

t-s>0 i->0 u 

which is the desired result. Finally suppose that ^o > and denote ke^^iO) = 
k{9o + 9)-k{eo). Trivially, 

ke^XO) = - (1 - e-^^)z/e„(rfx) with ue.idx) = e-^o^z/(dx) 
Jo 

and consider 

POO 

= e-'<^^Gix) - 00 J e-'-yG{y)dy (11) 

= fc,(^o) + G'(x) " (12) 



Line fllip is obtained by an integration by parts, whereas line f ll2p uses (jS]). 
We see easily from ( 1T^ that lim;j._j.oG6»o(2^)/log^ = ~^- Therefore we are in 
the same situation as in the proof of the first part with 60 = 0, and thus the 
proof is completed. ■ 

The following remark is useful to obtain a genralization of the three ex- 
amples presented in Section 3 for any i > 0. 

Remark 1 Suppose that {fXt)t>o is a family of infinitely divisible distribu- 
tions with fj,t* fJ's = fJ't+s, where t,s > 0. Assume that /ii fulfills the premises 
of Proposition 1 with G{x) ~ — £logx. Then obviously for any fixed t > 0, 
/ii also fulfills such premises with t£ replacing the role of L In all of the 
examples below, we have i = 1 and this remark shows how to get from them 
other examples with arbitrary i > 0. 

3 Examples 

Example 1 If fi{dx) = e~^l(^Q^oD){x)dx the corresponding infinitely divisible 
family is the gamma family with scale parameter 1. The Levy measure here 

l-(0,oo)( 



is u{dx) = e ""1 (0,00) (a;)— and i = 1. 



Example 2 A discrete example is 

00 



n=l 

[l/x] 



We have k{9) = ^^^^ ^(1 - e"^/"), G{x) = ^JJ/f 1/n ~ -logx z/a; ^ 
and i = l.The probability /i is the distribution of Yl'^=i^^ where X„ is 
Poisson distributed with mean 1/n and the {Xn)'^=i are independent. 

Example 3 Utilizing Example 2.2 in fJ^, consider the infinitely divisible dis- 
tribution fi on (0, 00) with Laplace transform defined on 6 > given by 
6 + 1 — \^(P + 26. Note that the densities of the corresponding EDM are 
Bessel densities related to a symmetric random walk (see Feller, 1971, pp. 
60-61). 

The related Levy measure of /i is 

dx 
v{dx)= iFi(l/2;l;-2x)l(o,oo)(a;) — , (13) 

Jb 

6 



where iFi{a; b; z) is the so called confluent entire function defined by Y^^=o nui) ^"' 
where (a)o = 1 and {a)^_^_^ = {a+n){a)n define the Pochhammer symbol (a)„. 
Therefore u has a density equivalent to 1/x when x — )■ which implies that 
G{x) ~ — logx. Proposition 1 is thus satisfied with i = 1. To check the 
correctness of (fT3l) observe that if k{9) = \og{9 + 1 — y/W+26) then 



f(x) = e 2"^ ^=1(0,00) (a;), ^(y) = ^=l(o,oo)(l/)- 



where 

Therefore the density of xv{dx), x > 0, is given by the convolution product 
/* (7, where by a change of variable y = tx and employing a Taylor expansion, 
one gets 

/ f{y)9{x-y)dy = - t- - dt= iFi(l/2;l;-2x). 



Let us fix t > 0. Recall (see [2]) that for 6* > the function [6 + 1-^/6'^ + 26) 
is the Laplace transform of the density ft{x) = -e~^ It{x)l(^o^oo){x) where th( 
Bessel function Itix) is 



It{x) = E 



^2n+i 



n\T(n + l + t) 2"+* ' 

n=0 ^ ' 

We fix now £ > 0. If l^i has density /« this implies that the density of 



U = Yf' is 



i2„,2t+l„ — V - ' 



It would be quite delicate to prove directly from this last formula that when 
t — )■ the law gt{u)du converges to the Pareto law 

idu 
as shown by our Proposition. 
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